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Abstract 

The calculation of the two-loop corrections to the three-jet production rate and to event shapes 
in electron-positron annihilation requires the computation of a number of two-loop four-point master 
integrals with one off-shell and three on-shell legs. Up to now, only those master integrals corresponding 
to planar topologies were known. In this paper, we compute the yet outstanding non-planar master 
integrals by solving differential equations in the external invariants which are fulfilled by these master 
integrals. We obtain the master integrals as expansions in e = (4 — d)/2, where d is the space- 
time dimension. The fully analytic results are expressed in terms of the two-dimensional harmonic 
polylogarithms already introduced in the evaluation of the planar topologies. 



1 Introduction 



Precision applications of particle physics phenomenology often demand theoretical predictions at the next- 
to-next-to-leading order in perturbation theory. Corrections to this order are known for many inclusive 
observables, such as total cross sections or sum rules. For 2^2 scattering and 1 — ► 3 decay processes, 
the calculation of next-to-next-to-leading order corrections is a yet outstanding task. One of the major 
ingredients for these calculations consists the two-loop virtual corrections to the corresponding four-point 
Feynman amplitudes. 

Using dimensional regularization [|]H with d = 4 — 2e dimensions as regulator for ultraviolet and 
infrared divergences, the large number of different integrals appearing in the two-loop Feynman amplitudes 
for 2 — > 2 scattering or 1 — > 3 decay processes can be reduced to a small number of master integrals. The 
techniques used in these reductions are integration-by-parts identities |^-^] and Lorentz invariance 0. A 
computer algorithm for the automatic reduction of all two- loop four-point integrals was described in |6j. 

For two-loop four-point functions with massless internal propagators and all legs on-shell, which are 
relevant for example to the next-to-next-to-leading order calculation of two-jet production at hadron collid- 
ers, all master integrals have been calculated over the past year [7|-|T2|. Very recently, these master integrals 
were already applied in the calculation of two-loop virtual corrections to Bhabha scattering in the limit 
of vanishing electron mass and to quark-quark scattering jl4|] . 

A different class of master integrals is required for the computation of next-to-next-to-leading order 
corrections to observables such as the three-jet production rate in electron-positron annihilation, two- 
plus-one-jet production in deep inelastic electron-proton scattering or vector-boson-plus-jet production at 
hadron colliders: two-loop four-point functions with massless internal propagators and one external leg 
off-shell. These functions involve one more scale than their on-shell counterparts, and are therefore more 
complicated and also more numerous. Up to now, only the master integrals for planar topologies [[l5|,[l6[ as 
well as one of the non- planar master integrals were known |^7| . In Jl^,|l^], Smirnov used a Mellin-Barnes 
countour integral technique to derive one planar and one non-planar master integral, obtaining analytic 
expressions for the divergent parts and an integral representation (one-dimensional in the planar case and 
three-dimensional in the non-planar case) for the finite part. Using the differential equation technique 
derived in we computed the full set of planar master integrals in Complementary work on the 

purely numerical evaluation of this type of master integrals has been presented recently by Binoth and 
Heinrich Q. 

Our results for the planar topologies [Hf ] were given in fully analytic form, in terms of a new class of 
functions, two-dimensional harmonic poly logarithms (2dHPL). The 2dHPL are an extension of the harmonic 
polylogarithms (HPL) of |nj. All HPL and 2dHPL that appear in the divergent parts of the planar 
master integrals have weight < 3 and can be related to the more commonly known Nielsen's generalized 
polylogarithms [^0[|2^] of suitable arguments. The functions of weight 4 appearing in the finite parts of the 
master integrals can all be represented, by the very definition, as one-dimensional integrals over Nielsen's 
polylogarithms of weight 3, hence of Nielsen's generalized polylogarithms of suitable arguments according 
to the above remark. A table with all relations is included in the appendix of ]l6[ ]. 

It is the purpose of the present paper to extend our work of towards the computation of non- 
planar four-point master integrals with one off-shell leg from their differential equations. In Section we 
discuss the use of the differential equation method in the computation of master integrals. We point out 
in particular the differences to the calculation of the planar master integrals. We then list our analytical 
results for the complete set of non-planar master integrals in Section ||. Our results are, like in the planar 
case, expressed in terms of 2dHPL. In Section || we also compare our results with the existing partial ones 
available in the literature. Finally, Section ^ contains conclusions and an outlook on future applications. 



2 Solving Differential Equations for Master Integrals 



The use of differential equations for the computation of master integrals was first suggested in |22| as a 
means of relating integrals with massive internal propagators to their massless counterparts. The method 
was developed in detail in (^3), where it was also extended towards differential equations in the external 
invariants. As a first application of this method, the two-loop sunrise diagram with arbitrary internal 
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masses was studied in [|24|. We have developed the differential equation formalism for two- loop four-point 
functions with massless internal propagators, three external legs on-shell and one external leg off-shell in [p| . 
We derived an algorithm for the automatic reduction, by means of computer algebra (using FORM |2q] 
and Maple |^6f), of any two- loop four-point integral to a small set of master integrals, for which we derived 
differential equations in the external invariants. 

Four-point functions depend on three linearly independent momenta: pi, P2 and p 3 . In our calculation, 
we take all these momenta on-shell (jpj — 0), while the fourth momentum p\2 3 = {p\ + P2 + P3) is taken 
off-shell. The kinematics of the four-point function is then fully described by specifying the values of the 
three Lorentz invariants s\ 2 = {pi +P2) 2 , S13 = (pi + P3) 2 and S23 = (P2 + p 3 ) 2 ■ 

Expressing the system of differential equations obtained in ^ for any master integral in the variables 
■S123 = S12 + S13 + S23, y — S13/S123 and z — S23/S123, we obtain a homogeneous equation in S123, and 
inhomogeneous equations in y and z. Since S123 is the only quantity carrying a mass dimension, the corre- 
sponding differential equation is nothing but the rescaling relation obtained by investigating the behaviour 
of the master integral under a rescaling of all external momenta by a constant factor. 

In the reduction of non-planar two-loop four-point functions with one off-shell leg, one obtains (in 
addition to master integrals corresponding to planar subtopologies) several genuinely non-planar master 
integrals, corresponding to three- and four-point topologies. The two master integrals corresponding to 
three-point topologies are the crossed vertex integral with one leg off-shell, two legs on-shell and the crossed 
vertex integral with two legs off-shell, one leg on-shell. The former obeys only a homogeneous differential 
equation, which is of no use for its complete computation, but it is otherwise relatively simple and already 
known in the literature J2?],[2£|]. We list its value in Section || for reference. For the latter vertex integral 
with two off-shell legs, one obtains one homogeneous and one inhomogeneous differential equation. The 
computation of this integral from its inhomogeneous differential equation is straightforward, following the 
procedure used for the planar integrals 16 1. 

Of the four non-planar four-point topologies with one off-shell leg, only one is reducible to planar 
subtopologies; the remaining three contain two master integrals each Q. For the six-propagator topology, 
we use the following basis of master integrals: 
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This choice of basis is not unambiguous. Our main motivation for this particular basis is the decoupling, 
in the d — > 4 limit, of the homogeneous part of the differential equations for each topology, which is crucial 
for the applicability of our algorithm to the solution of the differential equations outlined below. 

Quite in general, the analytic expression of all planar and non-planar amplitudes is found to consist 
of several terms, each equal to a prefactor (a simple rational polynomial in the external invariants) times 
a combination of 2dHPL multiplied by simple coefficients (integer numbers or rational coefficients). The 
major difference between the planar and non-planar master integrals is in the number of such terms. Indeed, 
all planar master integrals can be written as the product of a single prefactor times a Laurent series in 
e = (4 — d)/2 containing as coefficients only logarithms and polylogarithms of ratios of the invariants. It 
has become clear already in the on-shell case || that the non-planar master integrals do not take such a 
simple form. This expectation was confirmed for the off-shell non-planar master integrals in |Tg| ], where the 
1/e 4 and 1/e 3 terms of (§]|) and (|]| 

were calculated analytically. All these non-planar integrals contain 
the sum of several different prefactors, each multiplied by a different Laurent series in e. 

The algorithm employed in our calculation of the planar master integrals from their differential equa- 
tions ]l6f relied on the possibility of determining a unique rational prefactor for each master integral. Once 
this prefactor was found, the coefficients in the Laurent series were determined by inserting a general ansatz 
containing all possible harmonic polylogarithms, up to a given weight, and a subsequent matching of the 
individual terms. This procedure has to be modified to account for the more involved structure of the 
non- planar master integrals. 

With two master integrals T\(y, z, e) and T 2 (y, z, e) for a given topology, the differential equations in y 
take the following form, exact in e: 



d 
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T 2 (y,z,e) = R 21 (y,z,e)Ti(y,z,e) + R 22 (y,z,e)T 2 (y,z,e) + I 2 (y,z,e) 



(2.7) 



where the Rij(y, z, e) are rational functions of y, z and e, of order 1 in e; I\(y, z, e) and I 2 (y, z, e) are the 
inhomogeneous terms, containing the simpler and known master integrals of subtopologies of the topology 
under consideration. The first master integral T±(y, z,e) is always chosen to contain no scalar products 
in the numerator and all propagators only to unit power in the denominator, see (2.1),(2.3),( |2.5| ). The 
second master integral T 2 (y,z,e) can be chosen to contain one or more irreducible scalar products in the 
numerator, or one or more squared propagators in the denominator, or even a combination of both. The 
only requirement on T 2 (y,z,e) is its independence of T\(y,z,e) under integration-by-parts and Lorentz 
invariance identities. Any integral of the topology under consideration can then be expressed as a linear 
combination of T\ and T 2 plus simpler subtopologies. By exploiting the freedom to choose T 2 (y, z,e), we 
always succeeded in finding a T 2 {y 1 z, e) for which R 2 \(y, z, e) in (2.7) is of order e (this condition is fulfilled 
for the basis of master integrals ( |2.l| )-( |2~6| ) introduced above). As a consequence, the y-differential equation 
for T 2 (y, z, e), when Laurent-expanded in e, decouples from T\{y 1 z, e) and the coefficients of the expansions 
can be obtained from the differential equations by recursive quadrature. 



Let us illustrate the algorithm when Ti(y, z, e) and T 2 (y, z, e) are the amplitudes corresponding to (2.5) 
and (2.6), respectively. The explicit form of the equations, with d = 4 — 2e, in that case is 
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We now expand everything in Laurent series in e, up to terms of order 0, as 
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where the new inhomogeneous terms are 
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J%i{y,z) = -—T 14+1 (y,z) + —T 2 , i+ i(y,z) + l2,i{y,z); (2.11) 

and it is understood that quantities whose second index is greater than 4 vanish, as the greatest singularity 
of the Laurent expansion in e is 4, according to the analysis of the soft and collinear structure [^9|. The 
system is solved bottom up by starting from the coefficients of greatest singularity corresponding to i = 4, 
Tk,A{y, z), whose equations involve as inhomogeneous terms only the Ik,A{y, z), which are known. One then 
proceeds to the equations for Tk^(y,z), which involve as inhomogeneous terms Tk^{y,z) and Ik,a(y,z), 
which are known by now and so on until the T^fiiy, z) are obtained. 

According to the preceding general remarks, the equation for T 2 ,i{y, z) does not involve Ti i(y,z) (it 
involves T\ t i + i{y, z), which is known when considering the i-th term of the Laurent expansion); the homo- 
geneous part of the equation is almost trivial, having as a solution l/[y(l — y — z)]; by following the standard 
method of the variation of the constants, the solution of the inhomogeneous equation is immediately given 
by the quadrature formula 

T 2 Ay, z) - — r f Ay' y' (1 - y' - z)J 2 ,i(y', z) . (2.12) 

2/(1 - y - z) J 

The above result can be inserted into the equation for T\^{y, z), which is then solved in the same manner; 
starting from i = 4, the whole procedure can be iterated for a smaller value of i, till i = is reached (one 
could in principle arrive at any other desired value of i). Note that the whole process involves only the 
y-dependent denominators appearing in the differential equations in the variable y, namely 1/y, 1/(1 — y), 
l/(y + z) and 1/(1 — y — z). As the homogeneous terms are also equal to products of those denominators 
and 2dHPL functions, after careful partial fractioning of the y-dependent factors the integrations in ( 2.12| ) 



can always be evaluated analytically: when a square denominator, say 1/j/ 2 , occurs, they can be carried 
out almost trivially by integrating by parts, while when the integrand is a single denominator times any 
2dHPL function of some weight w, the primitive is nothing but another 2dHPL function of weight w + 1 
(see the Appendix of |ll| for more details on 2dHPL). 

The quadrature formula ( 2.12| ) is definite up to a constant of integration in y, so one still needs one 



constant per master integral per order in e. These constants of integration (which are functions of z) can be 
determined only by evaluating the master integrals for some boundary condition. In the planar case [fl6[ , 
we have used the fact that all planar master integrals, as well as their derivatives, are regular in the whole 
kincmatical plane, with the exception of two branch points at y — and z — 0. As a consequence, any 
factor (1 — y), (1 — y — z), or (y + z) appearing in the denominator of the homogeneous term of the differential 
equations for a master integral could be used to determine the boundary condition in y = 1, y = (1 — z), 
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or y = —z, respectively: multiplying the differential equation with one of these factors and taking the limit 
where the factor vanishes, the derivative and all the other terms not having that factor as a denominator 
drop out and one obtains a linear relation expressing the master integral in this special kinematical point in 
terms of the known master integrals of its subtopologies. The non-planar master integrals possess a different 
analytic structure: they have three branch points: in y = 0, z = and y = l — z. Consequently, only factors 
(1 — y) or (y + z) in the denominator of the homogeneous term can be used for the determination of the 
boundary conditions. It turns out that these boundary conditions alone are not sufficient for a complete 
determination of the constants of integration in all the master integrals: for each topology, one finds that 
one constant of integration per order in e remains unconstrained. We determine this (^-dependent) constant 
by solving the system of z-differential equations for each topology at the point y — 1 and requiring this 
solution to be regular at the point z = 1. Some care has to be taken in the analytic continuation of the 
master integrals across the cut in y = 1 — z to the point y = 1; details on the analytic continuation of the 
2dHPL can be found in the appendix of |l6| . 

Using the integration procedure described here, we have determined all non-planar two-loop four-point 
master integrals with one off-shell leg. The results are summarized in the following section. 



3 Master Integrals 



In the reduction of a non-planar two-loop four-point function with one off-shell leg, one will in general 
encounter master integrals with non-planar topology in combination with master integrals corresponding 
to planar subtopologies. These planar master integrals were all computed in jl6| . In this section, we tabulate 
all genuinely non-planar master integrals relevant to the computation of two-loop four-point functions with 
one off-shell leg. We classify the integrals, according to the number of different kinematical scales on which 
they depend, into one-scale integrals, two-scale integrals and three-scale integrals. 
The common normalization factor of all master integrals is 
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3.1 One-Scale Integrals 

The only non-planar one-scale integral is the crossed on-shell vertex integral, which was first computed 
in 1271 E8|. We list the result only for completeness: 
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3.2 Two-Scale Integrals 

The crossed vertex integral with two off-shell legs is the only non-planar two-scale master integral. This 
integral was, to the best of our knowledge, not known up to now. It fulfils one inhomogeneous differential 
equation in the ratio of the two external invariants, which can be employed for its computation. We obtain: 
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3.3 Three-Scale Integrals 

There are four different topologies of non-planar two-loop four-point functions with one off-shell leg, two 
topologies with six different propagators, and two others with seven different propagators. Only one of 
the six-propagator topologies is fully reducible to simpler (planar) subtopologies ||. All other topologies 
contain two master integrals each. The two master integrals for each topology fulfil a coupled system of 
inhomogeneous differential equations in y — S13/S123 and another coupled system in z = S23/S123. As 
explained in Section 0, solving one of these systems is already sufficient to compute the master integrals up 
to a constant of integration. The second system is used to constrain the boundary conditions and hence 
the constants of integration, as well as a check on the result. 

The choice of the basis of master integrals for each topology is described in Section ||. In particular, we 
always choose the integral without scalar products in the numerator and no squared propagator as one of 
the master integrals. This integral is symmetric under the interchange of two external momenta, which we 
take to be p\ and P2 ■ 

The two master integrals for the six-propagator topology are: 
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0, 0; z)//(l, 0; y) + 41/(0, 0; z)//(l - z, 0; y) + 4//(0, 0; z)H(l - z, 1 - z; y) 

0, 0, 0; y)//(0; z) - 4//(0, 0, 0; y)i/(l; z) - 41/(0, 0, 0; z)//(0; y) 

0, 0, 0; z)//(l - z; y) + 41/(0, 0, 0, 0; y) + 41/(0, 0, 0, 0; z) + 21/(0, 0, 0, 1; z) 

0, 0, 0, 1 - z; y) - 21/(0, 0, 1; z)H(0; y) - 61/(0, 0, 1; z)H(l; y) 

0, 0, 1; z)i/(l -z-y)- 21/(0, 0, 1, 0; y) - 61/(0, 0, 1, 0; z) + 41/(0, 0, 1, 1; z) 

0, 0, 1 - z; y)H(0; z) + 41/(0, 0, 1 - z; y)//(l; z) - 41/(0, 0, 1 - z, 0; y) 

0, 0, 1 - z, 1 - z; y) + 61/(0, 0, z; y)//(l; z) + 61/(0, 0, z, 1 - z; y) 



(3.13) 



0, 1; z)//(l - z, 0; y) - 41/(0, 1; z)//(l - z, 1 - z; y) + 81/(0, 1, 0; y)H(0; z) 

0, 1, 0; y)H(l; z) - 41/(0, 1, 0; z)H(0; y) - 101/(0, 1, 0; z)//(l; y) 

0, 1, 0; z)H(l -z;y)- 21/(0, 1, 0, 0; y) - 61/(0, 1, 0, 0; z) + 61/(0, 1, 0, 1; z) 

0, 1, 0, 1 - z; y) - 41/(0, 1, 1; z)//(0; y) - 41/(0, 1, 1; z)H{\ -z-y)- 91/(0, 1, 1, 0; y) 

0, 1, 1, 0; z) - 41/(0, 1, 1, 1; z) + 81/(0, 1, 1 - z, 0; y) - 61/(0, 1 - z; y)//(l, 0; z) 

0, 1 - z; y)H(l, 1; z) - 121/(0, 1 - z, 0; y)//(0; z) + 41/(0, 1 - z, 0; y)H(l; z) 

0, 1 - z, 0, 0; y) + 41/(0, 1 - z, 0, 1 - z; y) - 81/(0, 1 - z, 1, 0; y) 

0, 1 - z, 1 - z; y)H(0; z) - 41/(0, 1 - z, 1 - z; y)//(l; z) + 41/(0, 1 - z, 1 - z, 0; y) 

0. 1 - z, 1 - z, 1 - z; y) - 6//(0, 1 - z, z; y)H(l; z) - 61/(0, 1 - z, z, 1 - z; y) 
1; z)//(l - z, 0, 0; y) - 41/(1; z)H(l - z, 0, 1 - z; y) + 81/(1; z)H(l - z, 1, 0; y) 

1; z)//(l - z, 1 - z, 0; y) - 41/(1; z)H(l - z, 1 - z, 1 - z; y) - 81/(1, 0; y)H(l, 0; z) 

1, 0; y)l/(l, 1; z) - 101/(1, 0; z)i/(l, 1 - z; y) + 41/(1, 0; z)H(l - z, 0; y) 
l,0;z)H(l -z,l-z;y) - 81/(1, 0, 0; y)U(l; z) - 41/(1, 0, 0; z)//(l - z;y) 
1, 0, 0, 0;y) + 41/(1, 0,0,0;z) + 21/(1, 0,0, l;z) - 81/(1,0,0,1 - z;y) 

1, 0, 1; z)i/(l; y) + 41/(1, 0, 1; z)H(l - z; y) - 41/(1, 0, 1, 0; y) - 61/(1, 0, 1, 0; z) 

1, 0, 1, 1; z) + 81/(1, 0, 1 - z; y)H(l; z) - 81/(1, 0, 1 - z, 0; y) 

1, 0, 1 - z, 1 - z; y) + 61/(1, 0, z; y)//(l; z) + 61/(1, 0, z, 1 - z; y) 

1, 1; z)H(l - z, 0; y) - 41/(1, 1; z)//(l - z, 1 - z; y) - 101/(1, 1, 0; z)H(l; y) 
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with: 

fe.3A(y,z) 

f6.3,3{V> z ) 
h.3.2(Vi z ) 

fe.3,i{y,z) 



h.3,o(y,z) 



+61/(1, l,0;z)l/(l - z;y) + 6H(1, 1,0,0; y) - 61/(1, 1, 0, 0; z) + 61/(1, 1, 0, 1; z) 
-41/(1, 1, 1; z)H(l -z;y)- 31/(1, 1, 1, 0; y) - 51/(1, 1, 1, 0; z) - 41/(1, 1, 1, 1; z) 
+81/(1, 1 - z, 0; y)H(l;z) - 81/(1, 1 - z, 0, 0; y) + 81/(1, 1 - z, 0, 1 - z; y) 
-161/(1, 1 - z, 1,0; y) + 8H(l, 1 - z, 1 - z, 0; y) - 61/(1, 1 - z, z; y)H(l; z) 
-61/(1, 1 - z, z, 1 - z; y) - 41/(1 - z, 0, 0, 0; y) + 41/(1 - z, 0, 0, 1 - z; y) 
+2H(1 - z, 0, 1, 0; y) + 41/(1 - z, 0, 1 - z, 0; y) - 41/(1 - z, 0, 1 - z, 1 - z; y) 
-81/(1 - z, 1, 0, 0; y) + 8H(1 - z, 1, 0, 1 - z; y) - 61/(1 - z, 1, 1, 0; y) 
+81/(1 - z, 1, 1 - z, 0; y) + 41/(1 - z, 1 - z, 0, 0; y) - 41/(1 - z, 1 - z, 0, 1 - z; y) 
+81/(1 - z, 1 - z, 1, 0; y) - 41/(1 - z, 1 - z, 1 - z, 0; y) - 41/(1 - z, 1 - z, 1 - z, 1 - z; y) 
31tt 4 



+ - 



360 



+ c 3 



+ 101/(0; y) + 101/(0; z) - 131/(1; y) + 101/(1; z) + 2311(1 - z; y) 



7T 2 



41/(0; y)H(0; z) - 4H (0; y)H(l; z) - 101/ (0; z)H(l; y) + 611(0; z)H(l - z; y) 

-61/(0, 0; y) - 61/(0, 0; z) - 91/(0, 1; y) - 71/(0, 1; z) + 21/(0, 1 - z; y) - 101/(1; y)l/(l; z) 
+61/(1; z)H(l - z; y) - 121/(1, 0; y) - 61/(1, 0; z) - 31/(1, 1; y) - 71/(1, 1; z) 



-10// (1, 1 - z; y) + 6H(1 - z, 0; y) - 61/(1 - z, 1; y) + 121/(1 - z, 1 - z; y) 



Pl23-*~ 



P3- 



+-P2 



-pi 



\167T 2 / S12 + S13 



i=0 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



-1/(0; y)H(0; z) - 1/(0; z)H{\ - z; y) + 1/(0, 0; z) - 1/(0, 1; z) + — 

6 



(3.18) 

+21/(0; y)H(l, 0; z) + 21/(0; z)H(l - z, 0; y) - 2H(0; z)H(l - z, 1 - z; y) 

+21/(0, 0; y)l/(0; z) + 21/(0, 0; z)H(0; y) + 21/(0, 0; z)H(l - z;y) - 21/(0, 0, 0; z) 

+21/(0, 1, 0; y) + 31/(0, 1, 0; z) - 211(0, 1, 1; z) - 211(0, 1 - z; y)H(0; z) 

+211(0, z; y)H{\; z) + 21/(0, z, 1 - z; y) + 21/(1; z)H(l - z, z; y) + 21/(1, 0; z)H(l - z; y) 

+311(1, 1, 0; z) + 21/(1 - z, 1, 0; y) + 21/(1 - z, z, 1 - z; y) 



+ 



(3.19) 



+ 21/ (0; y) + 31/(0; z) + 31/(1; z) + 2//(l - z; y) 

-4//(0; y)H (1, 0, 0; z) - 4//(0; y)//(l, 1, 0; z) - 4//(0; z)H(l - z, 0, 0; y) 
-2//(0; z)//(l - z, 0, 1 - z; y) + 21/(0; z)H(l - z, 1, 0; y) - 61/(0; z)H(l - z, 1 - z, 0; y) 
-41/(0; z)//(l - z, 1 - z, 1 - z; y) - 41/(0, 0; y)//(0, 0; z) - 41/(0, 0; y)//(l, 0; z) 
+411(0, 0; z)H(0, 1 - z; y) - 411(0, 0; z)H{\ - z, 0; y) + 411(0, 0; z)H{\ - z, 1 - z; y) 
-411(0, 0, 0; y)H(0; z) - 411(0, 0, 0; z)H(0; y) - 41/(0, 0, 0; z)//(l - z; y) 
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+422 
-622 
-222 
-222 
-AH 
+222 
-AH 
-622 
+2H 
+AH 
-2H 
+AH 
-AH 
+2H 
-AH 
+AH 
-2H 
+AH 

120 

-AH 

-522 
-6H 



0, 0, 0, 0; z) + 2H(0, 0, 0, 1; z) - 2H(0, 0, 1; z)H(l -z;y)- 422(0, 0, 1, 0; y) 

0, 0, 1, 0; z) - 222(0, 0,1- z; y)H(0; z) + 2H(0, 0, z; y)H(l; z) + 2H(0, 0,z,l- z; y) 

0, 1; z)22(0, 1 - z; y) - 2H(0, 1; z)H(0, z; y) + 222(0, 1; z)H(l - z, 0; y) 

0, 1; z)H(l - z, z; y) + 21/(0, 1, 0; y)H(0; z) + 222(0, 1, 0; y)H(l; z) 

0, 1, 0; z)22(0; y) - 422(0, 1, 0; z)22(l - z; y) - 421(0, 1, 0, 0; y) - 622(0, 1, 0, 0; 2) 

0, 1, 0, 1; z) + 222(0, 1,0,1- z; y) - 2H(0, 1, 1, 0; y) - 522(0, 1, 1, 0; z) 

0, 1, 1, 1; z) + 222(0, 1,1-z, 0; y) - 622(0, 1 - z; y)H(l, 0; z) 

0, 1 - z, 0; y)22(0; z) - 422(0, 1 - z, 1, 0; y) - 422(0, 1 - z, 1 - z; y)2f (0; z) 

0, 1 - z, z; y)H{l; z) + 2H(0, 1 - z, z, 1 - z; y) + 422(0, z; y)22(l, 1; z) 

0, z, 1 - z; y)22(l; z) + AH(0, z, 1 - z, 1 - z; y) - 2H(0, z, z; y)H{l; z) 

0, z, z, 1 - z; y) + 422(1; z)22(l - z, 0, z; y) + 222(1; z)22(l - z, 1, 0; y) 

1; z)22(l - z, 1 - z, z; y) + 422(1; z)22(l - z, z, 1 - z; y) - 2H(1; z)H(l - z, z, z; y) 
1,0; z)H(l - z, 0; y) - 6H(1, 0; z)22(l - z, 1 - z; y) - 422(1, 0, 0; z)22(l - z; y) 

1, 0, 1; z)22(l - z; y) - 421(1, 0, 1, 0; z) + 4/2(1, 1; z)22(l - z, z; y) 

1, 1, 0; z)H(l -z-y)- 622(1, 1,0, 0; z) - 322(1, 1, 1, 0; z) - 622(1 - z, 0, 1, 0; y) 

1 - z, 0, z, 1 - z; y) - 422(1 - z, 1, 0, 0; y) + 222(1 - z, 1, 0, 1 - z; y) 

1 - z, 1, 1, 0; y) + 222(1 - z, 1, 1 - z, 0; y) - 622(1 - z, 1 - z, 1, 0; y) 

1 — z,l — z,z,l — z;y) + 422(1 — z,z,l — z,l — z;y) — 222(1 — z,z,z,l — z;y) 



+ c 3 



TV" 

IT 



- 422(0; y)22(0;z) 



1022(0; z) + 22(1; z) + 622(1 - z; y) 

0; y)22(l; z) - 422(0; z)H(l -z-y)- 422(0, 0; y) - 622(0, 0; z) - 222(0, 1; y) 
0, 1; z) - 422(0, 1 - z; y) - 622(1; z)H(l - z; y) - 422(1, 0; z) - 322(1, 1; z) 

1 - z, 0; y) - 222(1 - z, 1; y) - 622(1 - z, 1 - z; y) 



(3.20) 



The two master integrals for the first seven-propagator crossed topology can be expressed in the following 
form, which makes the pi «-» P2 interchange symmetry of (p. 3]) manifest: 



P123- 

P3- 





V 




A. 



Z>2 
Pi 



Vie^ 



(-S123 



r 2 '±r,{ 



S123 



^ e f ys? 2 si 3 s 2 3 
1 



S 12 S 13 



+ 



S? 2 fi 23 



/7.3,i 

/7.4,j 
/V.4,j 



gig S23 
Sl23 ' S123 

V s 123 ' S123/ 

( ^23 £13 
\Sl23 ' S123 



+ - 



l t ( S13 S23 \ 
/6.2,i , 



Si 2 Sl23 



\ s 123 Sl23 / 



[l - 4e + 16e 2 - 64e 3 ) 



(1 - 4e + 16e 2 ) + / 7 . 5)i 
(1 - 4e + 16e 2 ) + / 7 . 5)i 

I +0(e), 





S23 y 


\Sl23 


S123/. 




Sl3 \ 


\ s 123 


S123/ . 



(3.21) 



9 



and 



Pl23- 









A. 



P2 



,167T 2 y ^ \ S12S13S23 \ S 123 S123/ 



1 



S12S23 



\ s 123 ' s 123 / 



(8* 


Sl3 V 


\Sl23 


Sl23/_ 



+— /7.4,< f ^ (1 - 4e + 16e 2 ) ) + 0(e), (3.22) 

S12S13 \ s 123 5123/ / 



with 



fr. 3A (y,z) = 0, (3.23) 

fr.sAv,*) = -|, (3-24) 
h.z.Av^) = 3/7(0; y) + 3/7(0; z) (3.25) 
/r.3,i(tf.*0 = -6i/(0;y)i/(0;z)-6i/(0,0;y)-6i/(0,0;z)-6i/(l,0;y)-6i/(l,0;z)-7r 2 , (3.26) 
/ 7 .3,o(y, ^) = +12/7(0; y)/7(l, 0, z) + 12/7(0; z)H(l - z, 0, y) + 121/(0, 0; y)/7(0, z) 

+12/7(0, 0; z)H(0, y) + 12/7(0, 0, 0; y) + 121/(0, 0, 0; z) + 12/7(0, 1, 0; y) + 12//(0, 1, 0; z) 
+12//(1, 0; z)H{l -z;y) + 12/7(1, 0, 0; y) + 12/7(1, 0, 0; z) + 12/7(1, 1, 0; z) 
+12/7(1 - z, 1, 0; y) + 9( 3 + 2rr 2 [+77(0; y) + 77(0; z) + 77(1; z) + 7/(1 - z; y)] , (3.27) 
/ 7 . 4 , 4 (y,z) = 0, (3.28) 
/ 7 . 4 , 3 (y,z) = 0, (3.29) 
/ 7 .4, 2 (y,z) - -377(0; z)- 37/(1; z) - 37/(1 - z; y) , (3.30) 
/ 7 .4,i(y, z) = +677(0; y)/7(0; z) + 677(0; y)/7(l; z) + 677(0, 0; z) - 67/(0, 1; z) + 67/(0, 1 - z; y) 

-677(1; z)77(l - z; y) - 67/(1; z)//(z; y) + 67/(1, 0; z) - 67/(1, 1; z) + 67/(1 - z, 0; y) 
-677(1 - z, 1 - z; y) - 6/7(z, 1 - z; y) + rr 2 , (3.31) 
/ 7 .4,o(y, z) = -1277(0; y)/7(l, 0; z) + 1277(0; y)/7(l, 1; z) - 1277(0; z)77(l - z, 0; y) 

-1277(0, 0; y)/7(0; z) - 1277(0, 0; y)/7(l; z) - 1277(0, 0; z)/7(0; y) - 1277(0, 0, 0; z) 
-1277(0, 0, 1; z) - 1277(0, 0, 1 - z; y) + 127/(0, 1; z)77(0; y) - 127/(0, 1; z)77(l - z; y) 
-1277(0, 1, 0; z) - 1277(0, 1, 1; z) + 1277(0, 1 - z; y)/7(l; z) - 1277(0, 1 - z, 0; y) 
+1277(0, 1 - z, 1 - z; y) + 127/(0, z; y)/7(l; z) + 1277(0, z, 1 - z; y) 
+1277(1; z)77(l - z, 0; y) - 127/(1; z)//(l - z, 1 - z; y) - 127/(1; z)H(l - z, z; y) 
-1277(1; z)77(z, 1 - z; y) - 1277(1, 0; z)/Z(l - z; y) - 1277(1, 0, 0; z) - 127/(1, 0, 1; z) 
-1277(1, 1; z)/7(l - z; y) - 127/(1, 1; z)77(z; y) - 1277(1, 1,0; z) - 127/(1, 1, 1; z) 
-1277(1 - z, 0, 0; y) + 127/(1 - z, 0, 1 - z; y) - 2477(1 - z, 1, 0; y) + 127/(1 - z, 1 - z, 0; y) 
-1277(1 - z, 1 - z, 1 - z; y) - 1277(1 - z, z, 1 - z; y) - 1277(z, 1 - z, 1 - z; y) 
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-2tt 2 [+ff (0; y) + ff (0; z) + ff (1; z) + ff (1 - z; y)] , (3.32) 
f7.5,i{y,z) = -i , (3.33) 

/7.5,3(tf,*) = + iff(0;y)-^ff(0;z)-iff(l;z)-iff(l-z;y) , (3.34) 

/7.5, 2 (y, *) = +|ff(0; y)H(0; z) + iff(0; y)ff(l; z) - |ff(0; z)ff (1 - z; y) - iff(0, 0; y) + |ff(0, 0; z) 

+ |ff(0, 1; z) + iff (0, 1 - z; y) - iff (1; z)ff (1 - z; y) + 3ff (1; z)ff (z; y) - 2ff (1, 0; y) 

~H{1, 0; z) - iff (1, 1; z) + iff (1 - z, 0; y) - iff (1 - z, 1 - z; y) + 3ff (z, 1 - z;y) 

+y , (3-35) 
/7.6,i(t>. ^) = +3ff (0; y)ff (1, 0; z) + ff (0; y)ff (1, 1; z) - 4ff (0; z)ff (1, 0; y) + 3ff (0; z)ff (1 - z, 0; y) 
-3ff (0; z)ff (1 - z, 1 - z; y) - 6ff (0; z)ff (z, l-z;y)- 3ff (0, 0; y)ff (0; z) 
-ff (0, 0; y)ff (1; z) - 3ff (0, 0; z)ff (0; y) + 3ff (0, 0; z)ff (1 -z;y) + ff (0, 0, 0; y) 
-3ff (0, 0, 0; z) + 6ff (0, 0, 1; z) - ff (0, 0, 1 - z; y) - 3ff (0, 1; z)H (0; y) 
+3ff (0, 1; z)ff (1 - z; y) + 9ff (0, 1; z)ff (z; y) + ff (0, 1, 0; y) - 9ff (0, 1, 0; z) + 3ff (0, 1, 1; z) 
+3ff (0, 1 - z; y)ff (0; z) + ff (0, 1 - z; y)ff (1; z) - ff (0, 1 - z, 0; y) + ff (0, 1 - z, 1 - z; y) 
-6ff (0, z; y)ff (1; z) - 6ff (0, z, 1 - z; y) + ff (1; z)ff (1 - z, 0; y) 
-ff (1; z)H (1 - z, 1 - z; y) + 6ff (1; z)ff (1 - z, z; y) - 6ff (1; z)H (z, 0; y) 
+6ff (1; z)ff (z, 1 - z; y) + 15ff (1; z)ff (z, z; y) - 4ff (1, 0; y)H (1; z) 
-3ff (1, 0; z)ff (1 - z; y) - 6ff (1, 0; z)ff (z; y) + 4ff (1, 0, 0; y) + 3ff (1, 0, 0; z) 
+3ff (1, 0, 1; z) - 4ff (1, 0, 1 - z; y) - ff (1, 1; z)ff (1 - z; y) + 6ff (1, 1; z)ff (z; y) 
+3ff (1, 1, 0; y) - 3ff (1, 1, 0; z) - ff (1, 1, 1; z) - 4ff (1, 1 - z, 0; y) - ff (1 - z, 0, 0; y) 
+ff(l - z, 0, 1 - z; y) - 4ff(l - z, 1, 0; y) + ff (1 - z, 1 - z, 0; y) - ff (1 - z, 1 - z, 1 - z; y) 
+6ff (1 - z, z, 1 - z; y) - 6ff (z, 0, 1 - z; y) - 6ff (z, 1 - z, 0; y) + 6ff (z, 1 - z, 1 - z; y) 

+15ff (z, z,l-z;y) + iic3 + y[- 4ff (0; y) - 3ff (0; z) + 3ff (1; y) + 4ff (1; z) 

+4ff(l -z;y)] , (3.36) 
/7.5,o(y, ^) = -6ff (0; y)ff (1, 0, 0; z) - 6ff (0; y)H (1, 0, 1; z) + 6ff (0; y)ff (1, 1, 0; z) 

+2ff (0; y)ff (1, 1, 1; z) + 8ff (0; z)ff (1, 0, 0; y) - 2ff (0; z)ff (1, 0, 1 - z; y) 

+8ff (0; z)ff (1, 1 - z, 0; y) - 6ff (0; z)H (1 - z, 0, 0; y) + 6ff (0; z)H (1 - z, 0, 1 - z; y) 

-8ff (0; z)ff (1 - z, 1, 0; y) + 6ff (0; z)ff (1 - z, 1 - z, 0; y) 

-6ff (0; z)ff (1 - z, 1 - z, 1 - z; y) - 12ff (0; z)ff (1 - z, z, 1 - z; y) 

-6ff (0; z)ff (z, 0, 1 - z; y) + 8ff (0; z)ff (z, 1 - z, 0; y) - 12ff (0; z)ff (z, 1 - z, 1 - z; y) 

-6ff (0; z)ff (z, z, 1 - z; y) + 6ff (0, 0; y)H (0, 0; z) + 6ff (0, 0; y)ff (0, 1; z) 

-6ff (0, 0; y)ff (1, 0; z) - 2ff (0, 0; y)ff (1, 1; z) - 6ff (0, 0; z)ff (0, 1 - z; y) 



11 



+8/7(0, 0; z)H(l, 0; y) - 6/7(0, 0; z)/7(l - z, 0; y) + 6/7(0, 0; z)H(l -z,l-z;y) 

+12/7(0, 0; z)H(z, l-z;y) + 6/7(0, 0, 0; y)H(0; z) + 2H(0, 0, 0; y)H(l; z) 

+6/7(0, 0, 0; z)H(0; y) - 6/7(0, 0, 0; z)H(l - z;y) - 2H(0, 0, 0, 0; y) 

+6/7(0, 0, 0, 0; z) + 15/7(0, 0, 0, 1; z) + 2H(0, 0, 0, 1 - z; y) - 12/7(0, 0, 1; z)H(0; y) 

+6/7 (0, 0, 1; z)H(l; y) + 12/7(0, 0, 1; z)H(l - z;y) + 9/7(0, 0, 1; z)/7(z; y) 

+4/7(0, 0, 1, 0; y) - 6/7(0, 0, 1, 0; z) + 12/7(0, 0, 1, 1; z) - 6/7(0, 0, 1 - z; y)H(0; z) 

-2H(0, 0,1 -z; y)H{l;z) + 2H(0, 0,1 - z, 0; y) - 2H(0, 0, 1 - z, 1 - z; y) 

+12H(0, 0, z; y)H(l; z) + 12H(0, 0,z,l- z; y) - 6/7(0, 1; z)H(0, l-z;y) 

-12H(0, 1; z)/7(0, z; y) - 8/7(0, 1; z)H{l, 0; y) + 6/7(0, 1; z)H(l, l-z;y) 

-6/7(0, 1; z)/7(l - z, 0; y) + 6/7(0, 1; z)/7(l - z, 1 - z; y) + 18/7(0, 1; z)//(l - z, z; y) 

-12/7(0, 1; z)H{z, 0; y) + 12/7(0, 1; z)H{z, l-z;y) + 3/7(0, 1; z)/7(z, z; y) 

+8/7(0, 1, 0; y)/7(0; z) + 8/7(0, 1, 0; y)H(l; z) + 18/7(0, 1, 0; z)/7(0; y) 

+10/7(0, 1, 0; z)/7(l; y) - 18/7(0, 1, 0; z)/7 (1 - z;y) - 8/7(0, 1, 0; z)/7(z; y) 

-2/7(0, 1, 0, 0; y) + 18/7(0, 1, 0, 0; z) + 12H(0, 1, 0, 1; z) + 8/7(0, 1, 0, 1 - z; y) 

-6/7(0, 1, 1; z)/7(0; y) + 6/7(0, 1, 1; z)/7(l - z; y) + 18/7(0, 1, 1; z)H(z; y) 

-3/7(0, 1, 1, 0; y) + 6/7(0, 1, 1, 0; z) + 6/7(0, 1, 1, 1; z) + 8/7(0, 1, 1 - z, 0; y) 

+6/7(0, 1 - z; y)/7(l, 0; z) + 2H(0, 1 - z; y)H(l, 1; z) - 6/7(0, 1 - z, 0; y)/7 (0; z) 

-2/7(0, 1 - z, 0; y)/7(l; z) + 2/7(0, 1 - z, 0, 0; y) - 2H(0, 1 - z, 0, 1 - z; y) 

+2/7(0, 1 - z, 1, 0; y) + 6/7(0, 1 - z, 1 - z; y)/7(0; z) + 2/7(0, 1 - z, 1 - z; z) 

-2/7(0, 1 - z, 1 - z, 0; y) + 2H(0, 1 - z,l - z,l - z;y) - 12/7(0, 1 - z, z; y)H(l; z) 

-12/7(0, 1 - z, z, 1 - z; y) + 12/7(0, z; y)/7(l, 0; z) - 12/7(0, z; y)/7(l, 1; z) 

+6/7(0, z, 0; y)/7(l; z) + 6/7(0, z, 0, 1 - z; y) + 12/7(0, z, 1 - z; y)/7(0; z) 

-12/7(0, z, 1 - z; y)/7(l; z) + 6/7(0, z, 1 - z, 0; y) - 12H(0, z, 1 - z, 1 - z; y) 

-24/7(0, z, z; y)/7(l; z) - 24/7(0, z, z, 1 - z; y) - 2H{1; z)H(l - z, 0, 0; y) 

+2/7(1; z)/7(l - z, 0, 1 - z; y) - 12/7(1; z)/7(l - z, 0, z; y) - 8/7(1; z)H{l - z, 1, 0; y) 

+2/7(1; z)H{l -z,l-z,0;y)- 2H(1; z)H{l - z, 1 - z, 1 - z; y) 

+12/7(1; z)/7(l - z, 1 - z, z; y) - 12/7(1; z)/7(l - z, z, 0; y) 

+12/7(1; z)H{l -z,z,l-z;y) + 30/7(1; z)/7(l - z, z, z; y) + 12/7(1; z)/7(z, 0, 0; y) 
-12/7(1; z)H(z, 0, 1 - z; y) - 6/7 (1; z)/7(z, 0, z; y) - 12/7(1; z)/7(z, 1 - z, 0; y) 
+12/7(1; z)H{z, 1 - z, 1 - z; y) + 24//(l; z)//(z, 1 - z, z; y) - 6//(l; z)//(z, z, 0; y) 
+30/7(1; z)//(z, z, 1 - z; y) + 9//(l; z)//(z, z, z; y) + 8/7(1, 0; y)/7(l, 0; z) 
-8/7(1, 0; y)/7(l, 1; z) + 10/7(1, 0; z)/7(l, 1 - z; y) + 6/7(1, 0; z)/7(l - z, 0; y) 
-6/7(1, 0; z)/7(l - z, 1 - z; y) - 12/7(1, 0; z)/7(l - z, z; y) + 8/7(1, 0; z)/7(z, 0; y) 
+2/7(1, 0; z)/7(z, 1 - z; y) - 6/7(1, 0; z)/7(z, z; y) + 8/7(1, 0, 0; y)/7(l; z) 
+6/7(1, 0, 0; z)/7(l - z; y) + 12/7(1, 0, 0; z)/7(z; y) - 8/7(1, 0, 0, 0; y) 
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-6H(1, 0, 0, 0; z) + 12H{1, 0, 0, 1; z) + 8H(1, 0, 0, 1 - z; y) + 6H(1, 0, 1; z)H(l; y) 

+6H{1, 0, 1; z)H{\ -z;y) + 12iT(l, 0, 1; z)H(z; y) + 4ff(l, 0, 1, 0; y) 

-18# (1, 0, 1, 0; z) + 6if(l, 0, 1, 1; z) - 8H(1, 0,1 - z; y)H(l; z) + 8H(1, 0, 1 - z, 0; y) 

-8H(l, 0, 1 - z, 1 - 2r; y) - 6H(1, 0, z; y)H{l; z) - 6H{1, 0, z, 1 - z; y) 

+2iT(l, 1; z)if (1 - z, 0; y) - 2iJ(l, 1; z)H{\ - z, 1 - z; y) + 12iT(l, 1; z)iT(l - z, z; y) 

-12# (1, 1; z)ff(z, 0; y) + 12iT(l, 1; z)iJ(z, 1 - z; y) + 30H(1, 1; z)H(z, z; y) 

+10iJ(l, 1, 0; z)H(l; y) - 6H(1, 1, 0; z)H(l - z; y) + 2H(1, 1, 0; z)ff(z; y) 

-6H(1, 1, 0, 0; y) + 6H(l, 1, 0, 0; z) + 6ff(l, 1, 0, 1; z) - 2# (1, 1, 1; z)H(l - z; y) 

+12H (1, 1, 1; z)iJ(z; y) + 3^(1, 1, 1, 0; y) - 6JJ(1, 1, 1, 0; z) - 2JJ(1, 1, 1, 1; z) 

-8iJ(l, 1 - z, 0; y)H(l; z) + 8iJ (1, 1 - z, 0, 0; y) - 8H(1, 1 - z, 0, 1 - z; y) 

+16F(1, 1 - z, 1, 0; y) - 8#(1, 1 - z, 1 - z, 0; y) + 6H(1, 1 - z, z; y)if(l; z) 

+6iJ(l, 1 - z, z, 1 - z; y) + 2iJ(l - z, 0, 0, 0; y) - 2#(1 - z, 0, 0, 1 - z; y) 

+2iJ(l - z, 0, 1, 0; y) - 2H(1 - z, 0, 1 - z, 0; y) + 2#(1 - z, 0, 1 - z, 1 - z; y) 

-12# (1 - z, 0, z, 1 - z; y) + 8if(l - z, 1, 0, 0; y) - 8iJ(l - z, 1, 0, 1 - z; y) 

+6#(1 - z, 1, 1, 0; y) - 8iJ(l - z, 1, 1 - z, 0; y) - 2# (1 - z, 1 - z, 0, 0; y) 

+2H (1 - z, 1 - z, 0, 1 - z; y) - 8i?(l - z, 1 - z, 1, 0; y) + 2ff (1 - z, 1 - z, 1 - z, 0; y) 

-2iJ(l - z, 1 - z, 1 - z, 1 - z; y) + 12H (1 - z, 1 - z, z, 1 - z; y) 

-12H (1 - z, z, 0, 1 - z; y) - 12iJ(l - z, z, 1 - z, 0; y) + 12#(1 - z, z, 1 - z, 1 - z; y) 

+30iJ (1 - z, z, z, 1 - z; y) + 12iJ(z, 0, 0, 1 - z; y) - 14i/(z, 0, 1, 0; y) 

+12iJ(z, 0, 1 - z, 0; y) - 12#(z, 0, 1 - z, 1 - z; y) - 6iJ(z, 0, z, 1 - z; y) 

+12# (z, 1 - z, 0, 0; y) - 12iJ(z, 1 - z, 0, 1 - z; y) + UH(z, 1 - z, 1, 0; y) 

-12F(z, 1 - z, 1 - z, 0; y) + 12iJ(z, 1 - z, 1 - z, 1 - z; y) + 24iJ(z, 1 - z, z, 1 - z; y) 

-6#(z, z, 0, 1 - z; y) - 6ff (z, z, 1 - z, 0; y) + 30H (z, z, 1 - z, 1 - z; y) 



+9iJ(z, z, z, 1 - z; y) 



717T 4 C3 

240 



^ [ - 17H (0; y) + 9#(0; z) + 26if(l; y) + 17tf (1; z) 



+17if(l - z; y)] + — \ + 6H(0: y)H(0; z) - 8H(0; y)H(l; z) + 10.ff(0; z)H(l\ y) 
6 

-6iJ (0; z)H(l - z; y) + 8H(0, 0; y) + QH (0, 0; z) - 3iJ(0, 1; y) + 18if (0, 1; z) 
-8iT(0, 1 - z; y) + lOif (1; y)H(l; z) + 8H(1; z)H{\ - z; y) + 20#(1; z)fT (z; y) 
+12if(l, 0; y) - 67/(1, 0; z) + 3F(1, 1; y) + 8H(l, 1; z) + 10if(l, 1 - z; y) 
-8iT(l - z, 0; y) + 6H(1 - z, 1; y) + 8iJ(l - z, 1 - z; y) + 20ff (z, 1 - z; y)l . 



(3.37) 



The two master integrals for the second seven-propagator crossed topology are expressed as follows. 
Again, the p\ <-> P2 interchange symmetry of (|2.5|) becomes apparent: 



P123 
P3 



.X 









-P2 
-Pi 
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and 



16tt 2 / 



S12S13S123 \ s 123 «123/ 

S12S23S123 \ s 123 5123/ 

_! 1 ^ g / jSi3_ s 2 3 

Sl3(«13 + «23)si23 ' \ «123 ' §123 



1 



S23(S13 + S23)S123 '' \ s 123'si23 



h.J — ^—^+0{e), (3.38) 



P123 

P3 



■ WO) 






-*— 



-P2 



I671- 



^ £ YS12S13 \ s 123 S123/ 



— / 6 .2,i f — , ^) J + 0(e), (3.39) 

S13S123 \ s 123 5123/ J 

with 

/7.6,4(»,«) = -1, (3.40) 
/ 7 . 6 ,3(y,*) = 2ff(0;y) + 2H(0;z) , (3.41) 
/7.6,2(», «) = -4^(0; y)tf (0; z) - 4ff (0, 0; y) - 4H(0, 0; z) + 3H(1; z)H(l - z; y) - 3tf (1, 0; y) 

7T 2 

-3JT(1, 0; z) + 3i?(l, 1; z) + 3JT(1 - z, 1 - z; y) - — , (3.42) 

/ 7 .6,i(y, z) = +6tf (0; y)H(l, 0; z) - 4ff (0; y)H(l, 1; z) + 4tf (0; z)H(l, 0; y) + 2H(0; z)H(l - z, 0; y) 
-4ff(0; z)ff(l - z, 1 - z; y) + 8i?(0, 0; y)ff (0; z) + 8i?(0, 0; z)iJ(0; y) + 8H(0, 0, 0; y) 
+8H(0, 0, 0; z) + 6H(0, 1; z)iJ(l - z; y) + 6iJ(0, 1, 0; y) + 6ff(0, 1, 0; z) 
-2iJ(0, 1, 1; z) - 4if(0, 1 - z; y)i?(l; z) - 4£f(0, 1 - z, 1 - z; y) - 4if(l; z)i?(l - z, 0; y) 
+9iJ (1; z)H(l - z, 1 - z; y) + 10i?(l; z)i?(l - z, z; y) + 2H(1; z)H{z, 1 - z; y) 
+4ff(l, 0; y)H(l; z) + 2H(1, 0; z)ff(l - z; y) + 6# (1, 0, 0; y) + 6H(1, 0, 0; z) 
+10#(1, 0, 1; z) + 4#(1, 0, l-z;y) + 9H(l, 1; z)i?(l - z; y) + 2#(1, 1; z)H{z; y) 
+H(1, 1, 0; y) + 7JT(1, 1, 0; z) + 9if (1, 1, 1; z) + 4ff(l, 1 - z, 0; y) 
-4iJ(l - z, 0, 1 - z; y) + 6£f (1 - z, 1, 0; y) - 4if(l - z, 1 - z, 0; y) 
+9JT(1 - z, 1 - z, 1 - z; y) + 10/f (1 - z, z, 1 - z; y) + 2ff (z, 1 - z, 1 - z; y) 

+ I6C3 + y [ + 4F(0; y) + 4tf (0; z) + (1; y) + 3H(l;z) + 2H(l - z; y)] , (3.43) 
/ 7 .6,o(y, z) = -121/(0; y)ff (1, 0, 0; z) - 161/(0; y)ff (1, 0, 1; z) - 14ff (0; y)# (1, 1, 0; z) 
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-125(0; y)H{\, 1, 1; z) - 8H(0; z)H(l, 0, 0; y) - 25(0; z)H(l, 0, 1 - z; y) 

-85(0; z)H(l, 1 - z, 0; y) - 45(0; z)5(l - z, 0, 0; y) - 145(0; z)H(l - z, 0, 1 - z; y) 

-45(0; z)5(l - z, 1, 0; y) - 25(0; z)H(l - z, 1 - z, 0; y) 

-125(0; z)H(l - z, 1 - z, 1 - z; y) - 165(0, 0; y)5(0, 0; z) - 125(0, 0; y)H(l, 0; 0) 

+85(0, 0; y)5(l, 1; z) - 85(0, 0; z)5(l, 0; y) - 45(0, 0; z)5(l - z, 0; y) 

+85(0, 0; z)H(l -z,l-z;y)- 165(0, 0, 0; y)5 (0; z) - 165(0, 0, 0; z)5(0; y) 

-165 (0, 0, 0, 0; 2/) - 165(0, 0, 0, 0; z) - 125(0, 0, 1; z)H(l; y) + 145(0, 0, 1; z)H(l - z; y) 

-125 (0, 0, 1, 0; y) - 125(0, 0, 1, 0; z) + 45(0, 0, 1, 1; z) + 85(0, 0, 1 - z; y)5(l; z) 

+85(0, 0, 1 - z, 1 - z; y) - 125(0, 1; z)5(0, 1 - z; y) + 105 (0, 1; z)5(l, 0; y) 

-125(0, 1; z)H(l, 1 - z; y) - 145(0, 1; z)5(l - z, 0; y) + 135(0, 1; z)5(l - z, 1 - z; y) 

+65(0, 1; z)5(l - z, z; y) + 45(0, 1; z)5(z, 1 - z; y) - 45(0, 1, 0; y)5(0; z) 

-45 (0, 1, 0; y)5(l; z) - 125(0, 1, 0; z)5 (0; y) - 65(0, 1, 0; z)5(l; y) 

+25(0, 1, 0; z)H(l - z;y) - 125(0, 1, 0, 0; y) - 125(0, 1, 0, 0; z) - 45(0, 1, 0, 1; z) 

-45(0, 1, 0, 1 - z; y) + 11/7(0, 1, 1; z)5(l - z; y) - 45(0, 1, 1; z)5(z; y) 

+25(0, 1, 1, 0; y) - 65(0, 1, 1, 0; z) - 65(0, 1, 1, 1; z) - 45(0, 1, 1 - z, 0; y) 

-125(0, 1 - z; y)5 (1, 0; z) - 125(0, 1 - z; y)5 (1, 1; z) - 85(0, 1 - z, 0; y)5(0; z) 

+85(0, 1 - z, 0; y)5(l; z) + 85(0, 1 - z, 0, 1 - z; y) - 125(0, 1 - z, 1, 0; y) 

-125(0, 1 - z, 1 - z; y)5(l; z) + 85(0, 1 - z, 1 - z, 0; y) - 125(0, 1 - z, 1 - z, 1 - z; y) 

-125 (0, 1 - z, z; y)5 (1; z) - 125(0, 1 - z, z, 1 - z; y) + 8H(1; z)H(l - z, 0, 0; y) 

-125(1; z)5(l - z, 0,1 - z;y) - 45(l;z)5(l - z,l,0;y) - 125(1; z)5(l - z, 1 - z,0;y) 

+215(1; z)5(l - z, 1 - z, 1 - z; y) + 255(1; z)H(l - z, 1 - z, z; y) 

+235(1; z)5(l - z, z, 1 - z; y) + 65(1; z)5(l - z, z, z; y) + 65(1; z)5(z, 1 - z, 1 - z; y) 
+45(1; z)5(z, 1 - z, z; y) - 45(1; z)5(z, z, 1 - z; y) - 85(1, 0; y)5(l, 0; z) 
+105(1, 0; y)H(l, 1; z) - 65(1, 0; z)5(l, 1 - z; y) - 6H(l, 0; z)5(l - z, 0; y) 
-45(1, 0; z)H(l - z, 1 - z; y) - 85(1, 0, 0; y)5(l; z) - 45(1, 0, 0; z)H(l - z; y) 
-125(1, 0,0,0;y)- 125(1, 0,0,0; z) + 185(1,0, 0,1; z) - 85(1, 0,0,1- z;y) 
-125(1, 0, 1; z)5(l; y) + 91/(1, 0, 1; z)H(l -z;y) + 45(1, 0, 1; z)5(z; y) 
+25(1, 0, 1, 0; y) - 45(1, 0, 1, 0; z) + 215(1, 0, 1, 1; z) + 105(1, 0, 1 - z; y)5(l; z) 
-85(1, 0, 1 - z, 0; y) + 105(1, 0, 1 - z, 1 - z; y) + 125(1, 0, z; y)5(l; z) 
+ 125(1, 0, z, 1 - z; y) - 125(1, 1; z)5(l - z, 0; y) + 215(1, 1; z)5(l - z, 1 - z; y) 
+235(1, 1; z)5(l - z, z; y) + 65(1, 1; z)5(z, 1 - z; y) - 45(1, 1; z)5(z, z; y) 
-65(1, 1, 0; z)5(l; y) - 65(1, 1, 0; z)5(l - z; y) - 25(1, 1, 0, 0; y) - 145(1, 1, 0, 0; z) 
+75(1, 1, 0, 1; z) + 215(1, 1, 1; z)5(l - z; y) + 65(1, 1, 1; z)5(z; y) - 5(1, 1, 1, 0; y) 
-155(1, 1, 1, 0; z) + 215(1, 1, 1, 1; z) + 105(1, 1 - z, 0; y)5(l; z) - 85(1, 1 - z, 0, 0; y) 
+ 105(1, 1 - z, 0, 1 - z; y) - 185(1, 1 - z, 1, 0; y) + 105(1, 1 - z, 1 - z, 0; y) 
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-12H(1, l-z,z; y)H(l; z) - 12if(l, 1 - z, z,l - z;y) + 8H(1 - z, 0, 0, 1 - z; y) 

-6if(l - z, 0, 1, 0; y) + 8H(l -z,0,l-z, 0; y) - 12H{1 - z, 0, 1 - z, 1 - z; y) 

-12H{1 - z, 1, 0, 0; y) - 4H(1 - z, 1, 0, 1 - z; y) + 2JT(1 - z, 1, 1, 0; y) 

-4£T(1 - z, 1, 1 - z, 0; y) + 8H(1 - z, 1 - z, 0, 0; y) - 12#(1 - z, 1 - z, 0, 1 - z; y) 

+8H(l - z, 1 - z, 1, 0; y) - 12if (1 - z, 1 - z, 1 - z, 0; y) 

+2Lff(l - z, 1 - z, 1 - z, 1 - 2; y) + 25iT(l - z, 1 - z, z, 1 - z; y) 

+23-ff(l — z, z, 1 — z, 1 — z; y) + 6-ff (1 — z, z, z, 1 — z; y) 

+6ff (z, 1 - z, 1 - z, 1 - z; y) + AH(z, 1 - z, z, 1 - z; y) - 4if(z, z, 1 - z, 1 - z; y) 
2tt 4 

+ — + Cs [ - 32#(0; y) - 32tf (0; z) - llif(l; y) - 15if (1; z) - 4H(l - z; y)] 

+ y [ - 8H (0; y)tf (0; z) - 6if (0; y)# (1; z) - 6tf (0; z)H(l;y) - 81/(0, 0; y) - 8fT(0, 0; z) 

+2H(0, 1; y) - 6JT(0, 1; z) - 8tf (0, 1 - z; y) - 6H(l; y)H(l;z) - 4tf (1; z)ff (1 - z; y) 

-6tf(l, 0;y) -6#(1,0; z) - ff(l,l;y) - ll#(l,l;z) - 6H(1, 1 - z;y) 

+2ff (1 - z, 1; y) - 6tf (1 - z, 1 - z; y)] , (3.44) 

h.7,i{y,z) = ~, (3.45) 

h.7,3(y,z) = + b -H{0-y)+ i -H(Q ] z)+ l -H(hz)+ l -H{l-z-y) , (3.46) 

/7.7, 2 (y; *) = -3^(0; y)tf (0; z) - ff(0; y)# (1; z) + H(0; z)H(l - z; y) - 5tf (0, 0; y) - 3tf (0, 0; z) 
+3ff (0, 1; z) - (0, 1 - z; y) + AH(1; z)H(l - z; y) + 2tf (1; z)if (z; y) - 3H(l, 0; y) 

+4tf (1, 1; z) - £T(1 - z, 0; y) + 4tf (1 - z, 1 - z; y) + 2H(z, 1 - z; y) - ^ , (3.47) 

6 

/ 7 .7,i(y; z) = -8if(0; y)F(l, 1; z) + 4tf (0; z)tf (1, 0; y) + 4iJ(0; z)ff(z, 1 - z; y) + 6H(0, 0; y)tf (0; z) 
+2#(0, 0; y)H(l;z) + 6H(0, 0; z)H(0; y) - 2H(0, 0; z)H{\ - z; y) + 10i/(0, 0, 0; y) 
+6ff (0, 0, 0; z) + 9H(0, 0, 1; z) + 2H(0, 0, 1 - z; y) - 6iJ(0, 1; z)H(0; y) 
+9F (0, 1; z)tf (1 - z; y) + 7i? (0, 1; z)tf (z; y) + 6H(0, 1, 0; y) + 9tf (0, 1, 0; z) 
+6tf(0, 1, 1; z) - 2H(0, 1 - z; y)ff(0; z) - 8ff(0, 1 - z; y)H{l;z) + 2iJ(0, 1 - z, 0; y) 
-8i?(0, 1 - z, 1 - z; y) - 4ff(0, z; y)H(hz) - 4iJ(0, z, 1 - z; y) - 8#(1; z)iJ(l - z, 0; y) 
+lLff(l; z)tf(l - z, 1 - z; y) + 9ff(l; z)iJ(l - z, z; y) - 4H(l; z)H(z, 0; y) 
+6H(hz)H(z, 1 - z; y) + 3i/(l; z)iJ(z, z; y) + 4ff(l, 0; y)i?(l; z) + 2fl"(l, 0; z)iJ(l - z; y) 
+4ff (1, 0; z)H{z- y) + 6iJ(l, 0, 0; y) + 9ff(l, 0, 1; z) + 4fl"(l, 0, 1 - z; y) 
+ llff(l, 1; z)ff(l - z; y) + 6JT(1, 1; z)H(z; y) + H(l, 1, 0; y) + HJffl, 1, 1; z) 
+4ff(l, 1 - z, 0; y) + 2H(l - z, 0, 0; y) - 8ff (1 - z, 0, 1 - z; y) + 6JT(1 - z, 1, 0; y) 
-8i? (1 - z, 1 - z, 0; y) + llif (1 - z, 1 - z, 1 - z; y) + 9i?(l - z, z, 1 - z; y) 
-4tf(z, 0, 1 - z; y) - 4#(z, 1 - z, 0; y) + 6H(z, 1 - z, 1 - z; y) + 3if (z, z, 1 - z; y) 



16 



+ y C 3 + y [ + 2^(0; y) + 3H (0; z) + H(l;y) - 2H(l;z)] , (3.48) 

/ 7 .7,o(y; z) = -18^(0; y)if (1, 0, 1; z) - 221/(0; y)ff (1, 1, 1; z) - 8H(0: z)H(l, 0, 0; y) 

-2H(0; z)H(l, 0, 1 - z; y) - 8H(0; z)H(l, 1 - z, 0; y) - 2H(0; z)H(l -z,0,l-z; y) 

-2H(0; z)H(l -z,l-z,l-z;y) + 8H(0; z)H{z, 0, 1 - z; y) - 2H(0; z)H{z, 1 - z, 0; y) 

+AH (0; z)ff(z, 1 - z, 1 - z; y) - 12ff (0, 0; y)tf(0, 0; z) + 12H{0, 0; y)ff(0, 1; z) 

+ 16H (0, 0; y)tf(l, 1; z) + 4ff(0, 0; z)H(0, l-z;y)- 8H{0, 0; z)H(l, 0; y) 

-8tf (0, 0; z)H(z, l-z;y)- 12ff (0, 0, 0; y)ff (0; z) - 4ff(0, 0, 0; y)ff (1; z) 

-12JT(0, 0, 0; z)H(0; y) + 4ff(0, 0, 0; z)ff(l - z; y) - 20tf (0, 0, 0, 0; y) 

-12i?(0, 0, 0, 0; z) + 21£T(0, 0, 0, 1; z) - AH (0, 0, 0, 1 - z; y) - 18i/(0, 0, 1; z)H(0; y) 

-12H (0, 0, 1; z)ff(l; y) + 13#(0, 0, 1; z)H{\ -z;y) + 5H{0, 0, 1; z)H{z; y) 

-12fT(0, 0, 1, 0; y) + 18ff (0, 0, 1, 1; z) + 4ff(0, 0, 1 - z; y)ff (0; z) 

+ im (0, 0, 1 - z; y)tf(l; z) - 4ff(0, 0, 1 - z, 0; y) + 16#(0, 0, 1 - z, 1 - z; y) 

+8tf(0, 0, z; y)ff(l; z) + 8£T(0, 0, z, 1 - z; y) - 18tf(0, 1; z)H(0, 1 - z; y) 

-14ff(0, 1; z)JT(0, z; y) + lOff (0, 1; z)H(l, 0; y) - 12tf (0, 1; z)H(l, l-z;y) 

-18£f(0, 1; z)tf (1 - z, 0; y) + 21£f(0, 1; z)i7 (1 - z, 1 - z; y) - 3H(0, 1; z)iJ(l - z, z; y) 

-8tf (0, 1; z)H(z, 0; y) + 4if (0, 1; z)H(z, 1 - z; y) - 3ff (0, 1; z)ff (z, z; y) 

-8tf (0, 1, 0; y)H(0; z) - 8H(0, 1, 0; y)tf (1; z) - 18tf (0, 1, 0; z)tf (0; y) 

-6H(0, 1, 0; z)ff (1; y) + 2ff (0, 1, 0; (1 - z; y) + 2JT(0, 1, 0; z)H(z; y) 

-12H{Q, 1, 0, 0; y) - 18iJ(0, 1, 0, 0; z) - 8£T(0, 1, 0, 1 - z; y) - 12ff(0, 1, 1; z)iJ(0; y) 

+19iJ"(0, 1, 1; z)H{l -z;y) + 6tf(0, 1, 1; z)H(z; y) - 2H(0, 1, 1, 0; y) - 21£f(0, 1, 1, 0; z) 

+12tf(0, 1, 1, 1; z) - 8H(0, 1,1- z, 0; y) - 4ff(0, 1 - z; y)#(l, 0; z) 

-22H{Q, 1 - z; y)ff (1, 1; z) + 16£f(0, 1 - z, 0; y)H(l; z) - 4#(0, 1 - z, 0, 0; y) 

+16£T(0, 1 - z, 0, 1 - z; y) - 12ff(0, 1 - z, 1, 0; y) - 22if(0, 1 - z, 1 - z; y)i? (1; z) 

+16.ff(0, 1 - z, 1 - z, 0; y) - 22if(0, 1 - z, 1 - z, 1 - z; y) - 18H(0, 1 - z, z; y)i? (1; z) 

-18£f(0, 1 - z, z, 1 - z; y) - 8i?(0, z; y)H{l, 0; z) - 12i2"(0, z; y)ff (1, 1; z) 

+8ff(0, z, 0; y)H(l;z) + 8H(0, z, 0, 1 - z; y) - 8tf(0, z, 1 - z; y)iJ(0; z) 

-12£f(0, z, 1 - z; y)H(l: z) + 8ff(0, z, 1 - z, 0; y) - 12tf(0, z, 1 - z, 1 - z; y) 

-6tf(0, z, z; y)tf (1; z) - 6tf (0, z, z, 1 - z; y) + 16i2"(l; z)ff(l - z, 0, 0; y) 

-22H(l;z)H(l - z, 0, 1 - z; y) - 16^(1^)^(1 - z, 0, z; y) 

-22^(1^)^(1 - z, 1 - z, 0; y) + 25^(1^)^(1 - z, 1 - z, 1 - z; y) 

+23H(l;z)H(l - z, 1 - z, z; y) + 2Lff(l; z)i?(l - z, z, 1 - z; y) 

-3# (1; z)i7(l - z, z, z; y) + 8H(1; z)H{z, 0, 0; y) - 4ff(l; z)i7(z, 0, 1 - z; y) 

-16H(l;z)H(z, 0, z; y) - AH(1; z)H(z, 1 - z, 0; y) + 14ff(l; z)H{z, 1 - z, 1 - z; y) 

+2if(l; z)ff (z, z, 1 - z; y) - 3H (1; z)i/(z, z, z; y) - 8JT(1, 0; y)H(l, 0; z) 
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+105(1, 0;y)5(l,l;z) - 65(1, 0; z)5(l, 1 - z;y) - 2H(1, 0; z)H(z, 0; y) 

-65(1, 0; z)5 (z, 1 - z; y) - &H(l, 0, 0; y)H(l; z) - 45(1, 0, 0; z)H(l - z; y) 

-85(1, 0, 0; z)H(z; y) - 125(1, 0, 0, 0; y) + 155(1, 0, 0, 1; z) 

-85(1, 0, 0, 1 - z; y) - 125(1, 0, 1; z)5(l; y) + 215(1, 0, 1; z)5(l - z; y) 

+45(1, 0, 1; z)H{z; y) + 25(1, 0, 1, 0; y) + 215(1, 0, 1, 1; z) 

+105(1,0, 1 - z; y)H(l;z) - 85(1, 0,1 -z, 0; y) + 10H (1, 0, 1 - z, 1 - z; y) 

+125(1, 0, z; y)5(l; z) + 12ff(l, 0, z, 1 - z; y) - 22ff(l, 1; z)i7(l - z, 0; y) 

+25£f(l, 1; z)H(l - z, 1 - z; y) + 215(1, 1; z)H{l - z, z; y) - 45(1, 1; z)H{z, 0; y) 

+ 145(1, 1; z)5(z, 1 - z; y) + 25(1, 1; z)5(z, z; y) - 65(1, 1, 0; z)5(l; y) 

+45(1, 1, 0; z)5(l - z; y) - 65(1, 1, 0; z)5(z; y) - 25(1, 1, 0, 0; y) 

+215(1, 1, 0, 1; z) + 255(1, 1, 1; z)5(l - z; y) + 145(1, 1, 1; z)5(z; y) 

-5(1, 1, 1, 0; y) + 255(1, 1, 1, 1; z) + 105(1, 1 - z, 0; y)5(l; z) 

-85(1, 1 - z, 0, 0; y) + 105(1, 1 - z, 0, 1 - z; y) - 185(1, 1 - z, 1,0; y) 

+ 105(1, 1 - z, 1 - z, 0; y) - 125(1, 1 - z, z; y)5(l; z) - 125(1, 1 - z, z, 1 - z; y) 

-45(1 - z, 0, 0, 0; y) + 165(1 - z, 0, 0, 1 - z; y) - 125(1 - z, 0, 1, 0; y) 

+ 165(1 - z, 0, 1 - z, 0; y) - 225(1 - z, 0, 1 - z, 1 - z; y) - 165(1 - z, 0, z, 1 - z; y) 

-125(1 - z, 1, 0, 0; y) + 65(1 - z, 1, 1, 0; y) + 165(1 - z, 1 - z, 0, 0; y) 

-225(1 - z, 1 - z, 0, 1 - z; y) + 185(1 - z, 1 - z, 1, 0; y) - 225(1 - z, 1 - z, 1 - z, 0; y) 

+255(1 - z, 1 - z, 1 - z, 1 - z; y) + 235(1 - z, 1 - z, z, 1 - z; y) 

+215(1 - z, z, 1 - z, 1 - z; y) - 35(1 - z, z, z, 1 - z; y) + 85(z, 0, 0, 1 - z; y) 

-65(z, 0, 1, 0; y) + 85(z, 0, 1 - z, 0; y) - 45(z, 0, 1 - z, 1 - z; y) 

-165(z, 0, z, 1 - z; y) + 85(z, 1 - z, 0, 0; y) - 45(z, 1 - z, 0, 1 - z; y) 

+65(z, 1 - z, 1, 0; y) - 45(z, 1 - z, 1 - z, 0; y) + 145(2, 1 - z, 1 - z, 1 - z; y) 

317T 4 

+25(z, z, 1 - z, 1 - z; y) - 35(z, z, z, 1 - z; y) + 

+C 3 [ - 375(0; y) - 245(0; z) - 115(1; y) - 175(1; z) - 95(1 - z; y)] 

+ y [ - 65(0; y)5(0; z) + 45(0; y)5(l; z) - 65(0; z)5(l; y) + 25(0; z)5(l - z; y) 

-45(0, 0; y) - 65(0, 0; z) - 25(0, 1; y) - 155(0, 1; z) - 65(1; y)5(l; z) 

-25(1; z)H(z; y) - 65(1, 0; y) - 5(1, 1; y) - 45(1, 1; z) - 65(1, 1 - z; y) 

+45(1 - z, 0; y) + 65(1 - z, 1; y) - 45(1 - z, 1 - z; y) - 25(z, 1 - z; y)] , (3.49) 

/7.8,4(|/,*) = 0, (3.50) 

/ 7 . 8 , 3 (y,z) = 0, (3.51) 
/ 7 .8,2(y,*) - 0, (3.52) 
/ 7 . 8> i(y, z) = -25(0; y)5(l, 1; z) + 25(0; z)5(l - z, 1 - z; y) + 25(0, 1; z)5(l - z; y) + 25(0, 1, 1; z) 
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-2H 
+2H 
-2H 
-8H 
+AH 
+8H 
-AH 
+AH 
-8H 
-AH 
+AH 
+6H 
-AH 
+AH 
+AH 
-6H 
-AH 
+AH 
+AH 
+AH 
+2H 
-AH 
-AH 
+AH 
+2H 
-6H 
+AH 
-6H 
-8H 

+AH 
-AH 



1 - z; y)H(l; z) - 2H(0, l-z,l-z;y)- 2H(l;z)H(l - z, 0; y) 



1, 0; z)H(l -z;y) + 2H(1,0, 1; z) + 2H(1, 1,0; z) - 2H(1 -z,0,l-z; y) 
l-z,l-z,0;y) , (3.53) 
0; y)H{l, 0, 1; z) - AH(0; y)H(l, 1, 0; z) - 61/(0; y)H(l, 1, 1; z) 
0; z)H(l - z, 1, 0; y) + 6if (0; z)ff(l - z, 1 - 2, 1 - 2; y) 

0; z)H(l - z, z, 1 - z; y) - AH(Q; z)H{z, 1 - z, 1 - z; y) + AH(0, 0; y)H(l, 1; z) 

0, 0; z)if(l - z, 1 - z;y) + 12#(0, 0, 1; z)H{\ - z;y) - AH{Q, 0, 1, 1; z) 

0,0,1- z; y)H(l; z) + 4ff (0, 0, 1 - z, 1 - z; y) - AH(0, 1; z)H(0, l-z;y) 

0, 1; z)H(l - z, 0; y) + 6H(0, 1; z)if(l - z, 1 - z; y) + 8H(0, 1; z)H(l - z, z; y) 

0, 1; z)H(z, l-z-y)- 4ff (0, 1, 0; y)H(0; z) - AH(0, 1, 0; y)H(l; z) 

0, 1, 0; z)H{\ - z; y) + AH(0, 1, 0, 1; z) - 4ff(0, 1,0,1- z; y) 

0, 1, 1; z)H(l -z-y)- AH(0, 1, 1; z)H(z; y) - AH(0, 1, 1, 0; y) + 6if (0, 1, 1, 1; z) 

0, 1, 1 - z, 0; 2/) + AH(0, 1 - z; 0; z) - 6iJ(0, 1 - z; y)ff(l, 1; z) 

0, 1 - z, 0; y)H(0; z) + AH{0, 1 - z, 0; y)H(l;z) + 4ff(0, 1 - z, 0, 1 - z; y) 

0, 1 - z, 1 - z; y)if (Q; z) - 6if (0, 1 - z, 1 - z; y)if(l; z) + 4ff(0, 1 - z, 1 - z, 0;y) 

0, 1 - z, 1 - z, 1 - z; y) - 8i?(0, 1 - z, z; y)H(l;z) - 8H(0, 1 - z, z, 1 - z; y) 

0, z; y)i?(l, 1; z) - 4ff(0, z, 1 - z; y)ff(l; z) - 4iJ(0, z, 1 - z, 1 - z; y) 
l;z)ff(l - z,0,0;y) - 6if (1; z)H(l - z, 0,1 - z;y) - 8H{\-z)H{\ - z,0,z;y) 
l;z)#(l - z, 1, 0; y) - 6#(1; z)if(l - z, 1 - z, 0; y) - 8fT(l; z)ff(l - z, z, 0; y) 
1; z)if(z, 0, 1 - z; y) + AH(1; z)H(z, 1 - z, 0; y) - 4iJ(l, 0; z)H(l - z, 0; y) 

1, 0; z)H(l -z,l-z;y) + 8H (1, 0; z)ff(l - z, z; y) - 4# (1, 0; z)H{z, 1 - z; y) 
1, 0, 0; z)H(l - z; y) + 12if(l, 0, 0, 1; z) + 2H(l, 0, 1; z)# (1 - z; y) 

1, 0, 1; z)H(z; y) + 4#(1, 0, 1, 0; z) + 6# (1, 0, 1, 1; z) - 6ff(l, 1; z)iJ(l - z, 0; y) 

1, 1; z)H(z, 0; y) - 6ff(l, 1, 0; z)H(l -z-y)- AH(l, 1, 0; z)H(z; y) - AH(1, 1, 0, 0; z) 

1, 1, 0, 1; z) - 6H(1, 1, 1, 0; z) + 4#(1 - z, 0, 0, 1 - z; y) + 4#(1 - z, 0, 1 - z, 0; y) 

1 - z, 0, 1 - z, 1 - z; y) - 8i?(l - z, 0, z, 1 - z; y) + 4H(1 - z, 1, 0, 1 - z; y) 

1 - z, 1, 1, 0; y) + AH(1 - z, 1, 1 - z, 0; y) + 4ff (1 - z, 1 - z, 0, 0; y) 

1 - z, 1 - z, 0, 1 - z; y) + AH{1 - z, 1 - z, 1, 0; y) - 6i? (1 - z, 1 - z, 1 - z, 0; y) 

1 - z, z, 0, 1 - z; y) - 8H(1 - z, z, 1 - z, 0; y) + 4iT(z, 0, 1 - z, 1 - z; y) 

z, 1 - z, 0, 1 - z; y) + 4if(z, 1 - z, 1 - z, 0; y) + — [ - AH{0, 1; y) + 4ff(0, 1 - z; y) 

6 L 



1; z)tf (1 -z-y)- AH(1, 1; z) + 4£f (1 - z, 1; y)] 



(3.54) 



Both seven-propagator master integrals without scalar products in the numerator ( 3.21 ), ( |3.38| ) were 
first considered by Binoth and Heinrich [ fl8[ . Using a dedicated algebraic procedure to separate the structure 
of infrared divergencies, the 1/e 4 and 1/e 3 terms of these integrals were obtained analytically in ]l8| , the 
remaining terms were computed numerically order by order in e. We agree with the results for the 1/e 4 
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Smirnov has considered the problem of evaluating one of the seven-propagator 



and 1/e 3 terms of 
More recently, 

maste r integrals ( 3.21 ). The result, quoted in [p.7 |, contains explicit expressions for all divergent parts of 
( 3.21 ), which we confirm. The finite part of ( 3.21 ), however, is expressed as a three-dimensional Mellin- 
Barnes integral, which we are unable to transform into two-dimensional harmonic polylogarithms. A direct 
comparison with the finite part of (l7]] is therefore not possible. Based on general considerations about the 
analytic structure of the results obtained in fr^| , Smirnov moreover argues that two-dimensional harmonic 
polylogarithms could be insufficient to describe the finite part of ( 3.21 ) and calls for the introduction of 
three-dimensional harmonic polylogarithms - but that is not the shown by our results. 



4 Conclusion 

Two-loop four-point functions with one off-shell leg are an important ingredient for the calculation of next- 
to-next-to-leading order corrections to three-jet production and related observables in electron-positron 
annihilation. By exploiting integration-by-parts |^-^| and Lorentz invariance j(| identities, one can express 
the loop integrals appearing in these functions as a linear combination of a small set of master integrals, 
which are scalar functions of the external invariants. Their determination was up to now a major obstacle 
to further progress in next-to- next-to-leading order calculations. 

These master integrals fulfil inhomogeneous differential equations in their external invariants, which 
can be used for their computation. In this approach, one solves the differential equations by starting 
from an appropriate ansatz, which is determined from the homogeneous part of the equation. Once this 
ansatz is found, the equations are expanded in e = (4 — d)/2 and then solved by repeated quadratures, 
yielding a general solution that still contains free constants of integration. These constants are determined 
from boundary conditions to the same differential equations. Using this method, we have computed the 
master integrals for planar two-loop four-point functions with one off-shell leg and presented it in a previous 
publication fujfl . In the present work, we have now derived the corresponding non-planar master integrals 
as well. The results are expressed as a Laurent series in e, with coefficients containing 2dHPL, which 



were introduced in |16| as generalization of harmonic polylogarithms 19 and of Nielsen's generalized 
polylogarithms [^o|. All 2dHPL up to weight 3, as those appearing in the divergent parts of the master 
integrals, can be expressed in terms of generalized polylogarithms of complicated arguments, see |l6| , while 
the 2dHPL of weight 4 appearing in the finite parts of the master integrals involve a one-dimensional 
integral over a combination of generalized polylogarithms of weight 3. 

The results presented in this paper complete the calculation of master integrals needed for the compu- 
tation of two-loop four-point functions with one off-shell leg. These functions are a crucial ingredient for 
the virtual next-to-next-to-leading order corrections to processes such as three-jet production in electron- 
positron annihilation, two-plus-one-jet production in deep inelastic scattering, and vector-boson-plus-jet 
production at hadron colliders. The results given in |]l6| and in the present work correspond to the kine- 
matical situation of three-jet production in electron-positron annihilation, where all master integrals are 
real (up to a complex normalization factor). Using the prescription for the analytic continuation of two- 
dimensional harmonic polylogarithms derived in the appendix of ]l6| , our results can be continued into any 
other kinematic region of interest. 

Equipped with the complete set of master integrals, it should soon become possible to compute the two- 
loop virtual corrections to the above-mentioned processes. However, it must be kept in mind that these 
corrections form only part of a full next-to-next-to-leading order calculation, which also has to include 



the one-loop corrections to processes with one soft or collinear real parton 30-32J], as well as tree-level 
processes with two soft or collinear partons [ |33"| . Only after summing all these contributions (and including 
terms from the renormalization of parton distributions for processes with partons in the initial state), do the 
divergent terms cancel among one another. The remaining finite terms have to be combined into a numerical 
programme implementing the experimental definition of jet observables and event-shape variables. A first 
calculation involving the above features was presented for the case of photon-plus-one-jet final states in 
electron-positron annihilation in f33f| , thus demonstrating the feasibility of this type of calculations. 

Note added: After this paper was submitted for publication, we have carried out a detailed comparison 
with the purely numerical results of p8|. We find full agreement with the values quoted in 113] for the 
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master integrals (3.21) and (3.38) at specific points of the kinematical invariants, as well as with values 
for the master integral (3.9) provided by the authors of |[8). We would like to thank Thomas Binoth and 
Gudrun Heinrich for their assistance with this comparison. 
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